Partially filled stripes in the two dimensional Hubbard model: statics and dynamics 
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The internal structure of stripes in the two dimensional Hubbard model is studied by going beyond 
the Hartree-Fock approximation. Partially filled stripes, consistent with experimental observations, 
are stabilized by quantum fluctuations, included through the Configuration Interaction method. 
Hopping of short regions of the stripes in the transverse direction is comparable to the bare hopping 
element. The integrated value of ng compares well with experimental results. 



By now it is -well established that charged stripes are 
formed in a significant doping range of cuprate oxides 
[0-^. The existence of these stripes -was predicted, on 
the basis of mean field calculations, in advance of its 
observation, both in the one j^-Q] and the three bands 
Hubbard model ||]. This is one of the scarce theoretical 
results in the field of high-Tc superconductivity -which 
-was confirmed after its prediction. Interest in these cal- 
culations decreased, as it was generally understood that 
the Hartree-Fock approximation -was unable to obtain 
the partially filled stripes observed experimentally. Un- 
related calculations found stripes -with different fillings 
in the t-J model and, using different techniques, in 
the Hubbard model although other numerical cal- 
culations sho-w conflicting results ijll],^. Alternatively, 
it has been argued that stripes in doped Mott antiferro- 
magnets arise from a tendency to-wards phase separation, 
frustrated by electrostatic interactions • Stripes in 

the Hubbard model have also been analyzed -within slave- 
boson techniques, -which use the Hartree-Fock solutions 
as input [ p^ . 

In the present -work we sho-w that the mean field cal- 
culations, initially used to demonstrate the existence of 
stripes, can be systematically improved in order to study 
the partially filled stripes observed experimentally. In 
addition, they provide significant insight into the inter- 
nal structure of the stripes and their fluctuations in the 
transverse direction. The Hartree-Fock method can be 
considered a quasiclassical approximation to the spin and 
charge degrees of freedom. Their lo-w amplitude quan- 
tum fluctuations can be incorporated by using the Ran- 
dom Phase Approximation In addition, one needs 
to consider quantum tunneling processes bet-ween degen- 
erate, or nearly degenerate, Hartree-Fock solutions, -when 
there are many. This is achieved -with the Configuration 
Interaction method (CI), -widely used in quantum chem- 
istry jl^ . The combination of the Hartree-Fock and 
CI methods gives reasonable results even -when applied 
to one dimensional systems . The CI method restores 
the symmetries broken by the Hartree-Fock approxima- 
tion, and provides information on the quantum dynamics 



of the static solutions obtained in mean field, -which can 
be broadly classified into spin polarons or stripes. 

Previous mean field studies have focused on filled 
stripes , which tend to be the solutions with the lowest 
energy per hole in this approximation, especially for the 
values of J7/t ~ 4 used in the initial studies 0-0 . There 
are, however, self consistent solutions which describe par- 
tially filled stripes. These solutions have energies per hole 
comparable to those of the filled stripes for J7/i ~ 8 — 20. 

Typical stripe solutions obtained in this work within 
the mean field approximation and periodic boundary con- 
ditions, are shown in Fig.Q|, where, for simplicity, we 
show results obtained in a 5 x 4 cluster, and U = 8 
(hereafter we shall use t as the unit of energy). The 
filled stripe, like the ones considered in previous mean 
field calculations, is shown in Fig.|^ c), while Fig.||l| a) 
and b) correspond to half filled stripes, not reported ear- 
lier. These textures correspond to self consistent Hartree 
Fock solutions, and can be found in clusters of any size. 
A generalization of the solution shown in Fig.[0 b) to a 
16 X 16 cluster with 1/8 holes per site (32 holes) is shown 
in Fig.|Q. We have also calculated the energy required to 
create a hole in these configurations. In order to subtract 
the energy of the antiferromagnetic background we have 
averaged the energy per site obtained on 4 x 4, 4 x 6 and 
6x6 clusters, the result is 0.463 per site. The so obtained 
energies per hole are reported in Table I. It is interesting 
to note that the standard solution for the filled stripe 
shown in Fig. ||l[ c) is the best UHF solution only for suf- 
ficiently small U, while at intermediate U (including the 
case of t/ = 8 discussed here), there are better solutions 
which are inhomogeneous along the stripe direction, with 
tHF — — 1.081i, slightly smaller than the energies of the 
half-filled stripes. We also show in Table I the energy 
gain due to hybridization with displaced solutions of the 
same type (e^j). 
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FIG. 1. Spin and charge textures for the three stripe solutions discussed in the text. Solutions a) and b) correspond to hal 
filled stripes, while solution c) is a filled stripe. 



Half filled stripe 
Fig.[l] a) 


Half filled stripe 
Fig.[l] b) 


Filled stripe 
Fig.[l] c) 
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Table I. Energies, per hole, of the stripe solutions shown 
in Fig.[l] 

The best energies per hole that we have obtained, after 
combining, within the CI scheme, the previous configura- 
tions and others not shown in Fig.|Q] are eci — —1.573t 
(half filled stripes) and eci — —l.Sdit (filled stripes). 

The results in Table I have been obtained for a single 
stripe. They have uncertainities due to the cluster size 
dependence and the substraction of the antiferromagnetic 
background. In general, the energy differences are small, 
of a few hundredths of t. The stabilization of half filled 
stripes seems to be a general result, at least for values of 
U/t in the range 6-20. The value U/t = 8 corresponds 
to a t-J model with J/t = 1/2. The energy per hole 
that we obtain compares well with calculations for this 
model |§. For U/t = 12 (J/t = 1/3), we obtain eh ~ 
— 1.71 ±0.05t. This value is also in reasonable agreement 
with other calculations [ pT[ |. 

We have also calculated the hole-hole correlation, 
Cij = {{1 — ni){l — Tij)) for the half- filled stripe in a 
5x4 cluster. The optimal CI wavefunction, discussed 
above, was used in the calculation. Cij shows a maxi- 
mum (0.021) when the holes are separated by a vector 
Tij — Yi — Vj — (0, 2), that is when they are two lattice 
constants apart along the stripe direction. This feature of 
the hole-hole correlation differentiates stripes from other 
two hole configurations investigated within the t — J ||2^] 
or the Hubbard models p3| , p^ . We note, however, that 
Cij is also rather large (0.011) for the hole-hole separa- 
tion at which the results reported in those studies show 
the maximum hole-hole correlation, namely, r.^ = (1, 1). 

Our method is well suited to analyze quantum fluctu- 
ations of stripes in the transverse direction, as we can 



calculate the off diagonal terms of the full hamiltonian 
between stripe solutions which are displaced by an arbi- 
trary lattice constant. We find that the matrix element 
which describes the shifting of a stripe segment of length 
I by one lattice unit decays as t^e"'/'", where ts ^ t and 

Note that the two main configurations involved, Fig.||] 
a) and b), can be viewed as a superposition of site cen- 
tered and bond centered one dimensional domain walls 
|l9| . These configurations lead to a very accurate de- 
scription of the dynamics of a hole in the one dimensional 
Hubbard model jl^. The fact that the maximum hole- 
hole correlation is found along the stripe (see discussion 
above) implies an additional similarity between the stripe 
solution and one dimensional Luttinger liquids. Thus, 
the stripes can be viewed as a one dimensional system at 
quarter filling. 
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FIG. 2. Spin and charge textures for a 16 x 16 cluster with 
32 holes (a;=l/8, where x is the number of holes per site ) 
and U/t = 8. 

Because of the commensuration with the lattice, we 
find a gap in the charge spectrum. We estimate the value 
of this gap, Ach, beyond mean field theory, by calculating 
Ach = E„^+i + En^^i - 2En^ , around the optimal fiUing. 
We obtain, for U/t = 8, A^h = 0.7 ± O.lt, which is lower 
than the Hartree-Fock gap. For U/t = 12, we find A^h = 
0.5±0.1i, which compares well to the gap obtained using 
dynamical mean field techniques [Q. 

Upon doping, we expect that each electron leads to two 
charge 1 /2 solitons . The Hartree-Fock method is well 
suited to visualize charge fractionalization of this type 
p5[ , as shown in Fig. |^ a)] . The fact that the stripe is em- 
bedded in a two dimensional background allows for an al- 
ternative solution, the step also shown in Fig.^ b)]. 
This solution has a slightly lower energy. Thus, we expect 
that deviations from half filling in the stripes will lead to 
the formation of steps. We can calculate the quantum 
dynamics of the steps along the direction parallel to the 
stripe. Within the CI method, a step like the one shown 
in Fig..|^ b)] gives rise to a band of delocalized states of 
width 0.59t{U/t = 8) and 0.37t{U/t = 12). This band- 
width measures directly the delocalization of the step 
along the stripe due to quantum fiuctuations. The en- 
ergy gain due to this hybridization, 0.29< {U/t — 8) and 
0.18t{U/t ~ 12), further stabilizes the step. 



FIG. 3. Addition of a single hole to a half filled 
stripe. Hartree-Fock spin and charge textures calculated for a 
5 X 22 lattice with periodic boundary conditions and 12 holes 
(U/t = 8). a) Two kink solution, b) Step solution. 

The formation of steps is hindered when the stripes are 
sufficiently close. Experiments suggest that the closest 
distance between stripes is four lattice units, which is 
realized at 1/8 doping Near this filling, holes in 

the stripes will give rise to charge fractionalization and 
Luttinger liquid behavior p8| , p9[ . 

We have also computed the value of ng, from the so- 
lutions described above. They are shown in fig.Q, cal- 
culated in a 16 X 16 cluster, and obtained by integrating 
in energies over the top half of the valence band and 
the stripe states induced within the gap. Note that the 
charge and spin distributions are homogeneous in real 
space, as the CI method restores full translational sym- 
metry. The strong anisotropy, in k space, of n^^ reflects 
the existence of the stripe (see also [^,^ ) , and is con- 
sistent with the experimental results reported in ||27|. 
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FIG. 4. Values of rijj in a 16 x 16 cluster with periodic 
boundary conditions. 

The usefulness of the method used here for the study 
of half filled stripes suggests that a detailed description of 
the spin waves in the antiferromagnetic background is not 
essential Q. On the other hand, the scheme used here 
shows that there is a delicate balance between spin po- 
larons and stripes, which seems also to b e found in 

Montecarlo calculations of the t-J model \ 
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A detailed comparison of the relative stability of these 
phases cannot be done within the present method, and 
seems also difficult to achieve by other techniques. It 
would be interesting to analyze further the experimen- 
tal implications of this near degeneracy of qualitatively 
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different solutions. Finally, the present method suggests 
that inhomogeneous spin and charge textures, like stripes 
or spin polarons, can be viewed as microsocopic manifes- 
tations of phase separation, as the natural homogeneous 
solution which can be obtained in mean field has negative 
compressibility Q. 

In conclusion, we show that the half filled stripes ob- 
served in experiments in cuprates can be derived using 
the same methods which lead to the prediction of the 
existence of (filled) stripes in the Hubbard model [^-Q. 
The method allows us to describe the quantum trans- 
verse fluctuations of the stripes, which are significant 
(comparable to t) for segments of length of 4-6 lattice 
units. Finally, we present the distribution of ng, with 
good agreement with experimental results. 
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